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Abstract 
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From a data sample of 29058 r^ — > vr^Tr'^vr^Tr'^^'T- decays observed in the 
^ \ CLEO detector we derive a 95% confidence upper limit on the tau neutrino 

mass of 28 MeV. 
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I. MOTIVATION 

In the Standard Model the tau lepton and the tau neutrino form the third generation 
weak doublet of leptons. Most observations are consistent with zero mass for each of the 
three types of neutrino, and with the conservation of lepton number for each of e, fi, t 
separately. These suppositions, however, should be tested, especially in the light of recent 
indications [0,0 of oscillations among the neutrino species. 

There are model dependent limits on the possible values of the tau neutrino mass. The 
SuperKamiokande experiment measures the ratio of rates for i/^ and Ue from decay prod- 
ucts of particles produced in cosmic ray collisions in the atmosphere. The deficit in the u^ 
rate, along with its dependence upon neutrino zenith angle and energy, can be interpreted 
as due to oscillation of the u^ to z/^- with a maximal amplitude (sin 26' ~ 1) and a frequency 
determined by 5 x 10~^ < Am^ < 6 x 10~^ eV^ at 90% confidence. This would imply, if the 
Ur mass were much larger than the u^ mass, that 

0.02 < m^^ < 0.08 eV (90% c.l.). 

Astrophysical observations and cosmological theory limit the energy density of the uni- 
verse, thus restricting the sum of stable neutrino masses [^. This leads to the limit 

m^^ < 24 eV. 

This limit not only depends on z/,- being stable, but also varies with the value of the Hubble 
constant and other inputs. 

The possible effects of the tau neutrino on big-bang nucleosynthesis lead to either a low 
mass Ut- or to a decaying u^- at higher mass [Q: 

m^^ < 0.37 MeV, or m^^ > 18 MeV. 

The width of the excluded region increases with the lifetime of the z/^ and also depends on 
the assumed abundance of light elements. 

The Standard Model relations among the r and n masses and lifetimes and some r 
branching fractions ^ imply 

m^^ < 48 MeV. 

The popular see-saw mechanism for generating neutrino masses postulates the relation 

m^^ : m^^ : m^^ = ml : mf^ : ml. 
This would imply rather weak limits on the Pr mass: 

rrii,^ < 180 MeV, from rni^,^ < 15 eV, 

m^^ < 48 MeV, from m^^^ < 0.17 MeV. 

The model dependences in all of these limits on the z/,- mass argue for a more direct 
measurement. The observation of the decay of accelerator produced taus along with mea- 
surements of the energy and momentum of the detectable daughter products (all but the 



Ut-) can constrain the possible values for rriu^, especially in cases when the effective mass of 
the detected particles is close to m,-. Because the cross section for t^t~ production is up 
to 15% of the total e~^e~ annihilation cross section, an electron-positron collider is a natural 
choice for the source of taus. 

One looks for a hadronic decay mode with only one neutrino in the final state. Most 
of the previous measurements (see Table |) have been made with the higher multiplicity 
decay modes in which the effective mass of the hadrons is more likely to be close to the 
kinematic limit m,- with the maximum sensitivity to m,^^. All such decay modes are strongly 
phase-space suppressed, however, and the branching fractions are very low. An alternative 
strategy, which we use in the present measurement, is to pick a decay mode with a lower 
hadron multiplicity but with a much larger branching fraction. Although the four-pion decay, 
r^ -^ 7r^7r+7r^7r°z/T-, relative to the five-pion decay, produces a smaller proportion of events in 
which the effective hadronic mass is close to rrir, the branching fraction is 4.2%, as compared 
to ~ 0.1% for the higher multiplicity modes. 

II. DATA SAMPLE AND EVENT SELECTION 

The experiment was performed using the Cornell Electron-positron Storage Ring and the 
CLEO II detector, described elsewhere ^. Charged particle tracks were reconstructed in 
three nested cylindrical drift chambers in a solenoid field of 1.5 Tesla. The mean-squared 
resolution in momentum transverse to the beam was {5pt/ptY — 0.005^ + (0.0015 GeV~^ x 
PtY- Photon and electron showers were detected over 98% of 47r steradians in an array of 
7800 Csl scintillation counters with an energy resolution of 6E/E = 0.0035/ii^°'^^ + 0.019 — 
O.OOIE' {E in GeV) in the central region of polar angle, 45° < 9 < 135°. Ionization, time of 
flight, and shower energy aided in lepton identification. 

The data used in the present analysis were from 4.75 fb^^ of accumulated luminosity, 
two thirds at 10.58 GeV and one third at 10.52 GeV e^e~ center of mass energies. This 
corresponds to 4.3 million t^t^ pairs produced. 

We determined event selection criteria using Monte Carlo simulated signal and back- 
ground data samples. We selected the one-versus-three track topology with zero net charge, 
that is, events containing a three-charged-track tau decay candidate tagged by a single- 
prong decay in the opposite hemisphere. Tracks were accepted in the polar angle range 
I cos6'| < 0.9. The three signal tracks each had to have pr > 0.019£^feeam and had to fail elec- 
tron identification criteria [§. The tag track had to have momentum greater than 0.047ii^feeam 
and had to be consistent with one of four possible decay modes: 

or the charge conjugates. This resulted in a data sample of 813 thousand tagged three-prong 
tau decay candidates. 



On the three-prong side each event was required to have a 7r°, defined as two Csl 
calorimeter showers in the polar angle range | cos6'| < 0.71, not matched to charged tracks, 
with lateral shower profiles consistent with photons, and with effective mass in the range 
120 < m^^ < 150 MeV. The 7r° in the decay of the p"*" in the r"*" — >■ p'^T^ tag channel also had 
to satisfy these requirements, along with E^ > 50 MeV for each photon and ?7i(7r+7r°) < 1.5 
GeV. Non-photon calorimeter showers can make false vr*^ candidates. These typically origi- 
nate from nuclear interactions of the charged pions in the Csl crystals producing secondaries 
isolated from any charged track. If several vr" candidates were found on the three-prong side 
of an event, we kept only the one with the highest energy. We rejected events with an extra 
shower of more than 300 MeV or, if photon-like in lateral shower shape, of more than 100 
MeV. The 7r° selection cuts reduced the data sample to 31,305 events. 

In order to minimize background from the two-photon process, for example e~^e~ -^ 
e"'"e~7*7*, 7*7* — > 7r"'"7r"'"7r~7r~7r'^ with the final e"*" and e~ escaping detection at small angles 
to the beam, we rejected events in which the net event momentum transverse to the beam 
was less than 150 MeV. The final data set included 29,058 events. 

No particle identification information was used on the charged particles. That is, decay 
modes in which a K^ substitutes for a vr^ were considered part of the r~ -^ h~h'^h~Ti^Vj. 
signal. These modes contributed about 5% of the signal (see Table ||). The four-pion final 
state includes K^n'Ti^Vr, K^ -^ tt+tt", at 3% of the total. 

The bulk of the signal is r^ -^ -K^Ti^'K^'K^Vr (and charge conjugate). In about 53% of 
the events there is a tt+tt^tt'^ combination (Fig. |l]) with a mass consistent with the uo (there 
are two possibilities per event). The two-pion mass spectra (Fig. |^) for events with no uj 
show p peaks in 7r"'"7r° (21% per event), 7r~7r° (17% per event), and vr+vr" (< 2% per event). 
The overall four-pion mass spectrum (Fig. |^) has a broad maximum around 1.2 to 1.4 GeV 
for the ujTT~ events, and a peak at 1.4 GeV for the rest of the events. There is no obvious 
resonance structure in m4^, although the two four-pion spectra fit well each to a sum of 
p(770), p(1450), and p(1700) with adjustable relative amplitudes and phases. 

III. ANALYSIS 

Conservation of energy and momentum imply that 



= mi + mjj- 2EbeamEH + 2V^L,n " ^l\J Eji - mi COS Oht- 

The beam energy Ef,eam and the energy Eh and effective mass rriH of the hadronic (four- 
pion) system are measured in each event, and m^ = 1777.05lQ;|g MeV is known P,p!0|. Thus 
if we fix rriiy^, for each allowed value of the scaled hadronic mass x = mH/TriT- there is a 
range of kinematically allowed values of the scaled hadronic energy y = En/Ebeam, where 
the limits are obtained by taking the last term to be at its cos 6ht = =tl limits. Figure ^ 
shows the distribution of the data in x,y and the boundary curves for two values of m^^. 
Even though there is background outside the kinematically allowed region, it is clear that 
the two-dimensional distribution of the data is sensitive to the value of m,^^ . More precisely, 
the likelihood of the observed x, y event distribution, including background, plotted for 



various assumed rrii^,^ values, can give information on which m,^^ values are consistent with 
experiment. 

We define the likelihood for an individual event observed at Xi, t/i as the probability 
density P{xi,yi\m,y^) of observing such an event assuming m,y^ to be the neutrino mass. 
P{xi,yi\my^) contains terms for signal and backgrounds. The likelihood for the entire data 
sample is then the product of the event likelihoods: 

AT 

C{my^) = Y[P{xi,yi\m^^). 

Signal Likelihood. We first discuss the signal contribution to the single-event probability 
density P{xi,yi\m^^). It can be expressed as the product of the spectral function (decay 
probability density) JF and detection efficiency e at the true x, y, convolved with the experi- 
mental resolution function 7?., derived from data for that event: 



Psigixi, yi\m^r) = / ^{^1 y\0)w{x, y\m,y^)e{x, y)7li{xi -x,y^- y)dxdy. 

For convenience we have expressed J^(x, y\Tny^) in terms of the spectral function for my^ = 
and a weight function that takes account of the dependence of J-' on rriu^ . The weight function 
w is determined from the known effect of a non-zero neutrino mass on the phase space and 
the kinematic boundary. It is zero outside the allowed region. 

The spectral function JF(x, y|0) is obtained by adjusting a physics motivated 14-parameter 
function [§ to match distributions in the simulated data and the real data over the range 
X < 0.925, where we have verified by Monte Carlo that the choice of JF does not bias the 
determination of rriy^. The function so determined is then used in the entire x < 1 range. 
It includes adjustable amplitudes and phases for cuvr and for pvrvr in all charge combinations. 
The UTi and pirn mass spectra are each a superposition of p(770), p(1450), and p(1700) 
resonances. The masses and widths of the resonances are fixed. We adjust the parameters 
by comparing distributions in simulated data, including the effects of detector acceptance 
and resolution, and in real data with estimated tau and non-tau backgrounds subtracted. 
Figures and ^ illustrate the goodness of the fit. 

We compute the resolution function TZi separately for each event. The scale of the 
spreading in x — Xi and y — yi, including x,y correlation, is obtained by propagating the 
resolution error matrices from the individual track and shower fits. We obtain a parametrized 
non-Gaussian shape from Monte Carlo. The width and shape lead to distributions that match 
data for the reconstructed mass of vr'' -^ 77, Kg -^ tt+tt^, and D — > Ktt. The projected 
distributions of r.m.s. resolutions in hadronic mass and energy peak at 11 and 17 MeV in 
the region near the kinematic endpoint. By using the resolution function appropriate for 
each event instead of an averaged one, we diminish the effect of fluctuations from poorly 
measured events near the kinematic boundary. 



We evaluate the integral for each event by Monte Carlo, using a GEANT |]Tl| simulation of 
physical processes in the CLEO detector. That is, we flrst generate about 1 million simulated 
signal events using the KORALB event generator |T^, tagged r~ — »• 7r~7r"'"7r~7r°^',-, according 



to the distribution J-'{x, y\0). Then for each observed real event i and for each assumed value 



of rrii,^ we form the following sum over all the Monte Carlo events that contribute: 



^sig Kp^i ) Vi I T^u-r , 



Ej=i^'wixj,yj\m^,] 



We take account of the efficiency factor e by omitting the Monte Carlo events that are not 
detected by the simulated CLEO detector or recognized by the event selection criteria. The 
Monte Carlo integration technique enables us to include in P^ig the effect of initial state 
radiation, e+e~ -^ r^r^'y. Radiative events have a lower effective Ef,eam causing them to be 
produced with lower apparent y. Some of them can be seen in Fig. ^ below the lower no- 
radiation kinematic limit. The Monte Carlo also includes the appropriate number of events 
from the Ks — >■ 7r"'"7r~ and i^^-for-vr^ substitution modes (see Table 0). 

Background Likelihood. We distinguish three types of significant background: (a) events 
from the two-photon process that are not eliminated by our transverse momentum cut, (b) 
T^T~ events that do not contain our signal modes, and (c) non-r hadronic events from 
e'^e" -^ qq {q = u, d, s, c). 

The two-photon events, such as e^e^ -^ e^e^7*7*, 7*7* -^ Tr+vr+Tr^Tr^vr'^ in which the 
hadronic state has enough transverse momentum to be accepted, form a background at low 
scaled hadron energy y that is difficult to model reliably. Since this kinematic region is 
insensitive to m^^, the best strategy is to eliminate it from the likelihood fit. The detec- 
tor efficiency is also less accurately modeled at low x and low y, so we restrict the fit to 
x>0.7, y>0.7. Within this region the two-photon background can be neglected, and the 
detailed choice of boundary has no infiuence on the m^^ limit. The cut reduces the number 
of data events used in the fit to 16,577. 

Tau decays of higher or lower multiplicity can masquerade as our signal mode if particles 
escape undetected and/or secondaries in the Csl crystal array are misinterpreted as photons 
from a vr*^. We evaluate these and other mis-reconstruction effects by Monte Carlo simu- 
lation of the response of the CLEO detector to 12 million t~^t~ events generated with the 
known branching fractions. Of the accepted data events in the fit region, 7.3% are from tau 
background. They are mainly r^ -^ -k^-k^-k^Vj. with a spurious vr^. 

Although most hadronic e"'"e^ -^ qq annihilation events are rejected by our one-versus- 
three charged track criterion, some of them can survive. The qq contamination in our data 



sample was evaluated by a 36-million event Monte Carlo simulation using the LUND ||13 
generator. The simulation has been extensively tuned to produce results that agree with 
experiment. In particular, we have verified the agreement between data and Monte Carlo 
for the events satisfying the one-versus-three topology but having a tag with an energy that 
could not come from r decay, and for the events that have x values well above the kinematic 
limit for the signal r mode. The qq background accounts for 3.1% of the accepted data 
events in the fit region. 

The calculation of the background contribution P^fc^ to the individual event likelihood is 
similar to the calculation of Psig. However, since the expected background event distribution 
(the analog of JF x e) can be expressed only in terms of observed Xi,yi, it is not appropriate 
to integrate over the experimental resolution function; its effect is already contained in the 
distribution. As the distribution is not an analytic function, but a collection of simulated 
events with a rather smooth distribution, we approximate the value of Pbkg{xi,yi) for the 



i-th data event by the number (appropriately weighted) of t^t and qq background Monte 
Carlo events per unit area in the x, y vicinity of Xi,yi. This Pbkg is of course independent of 

We sum the signal and background likelihoods for each event and take the product over 
all events in the region x > 0.7, y > 0.7 to form the net likelihood C This is repeated for 
a range of assumed m^^ values to obtain £(mi^^).[] The overall normalization of C{mi,^) is 
arbitrary. 

IV. RESULTS 

Figure ^(a) shows the likelihood as a function of assumed neutrino mass. The integral 
under the curve beyond 22 MeV is 5% of the total. Before interpreting this as a 95% upper 
limit on m,^^, however, we have to consider systematic uncertainties that could affect the 
limit. 

The CLEO charged particle momentum measurement scale is uncertain by about 0.05%, 
and the vr*^ energy scale is uncertain by 0.25%. This shows up as a potential mismatch 
between Monte Carlo and data, resulting in a distortion in C{my^). Variations in scale 
of these magnitudes cause a shift of 5.0 MeV in the 95% limit, when the two effects are 
combined in quadrature. 

The four-pion spectral function JF was determined by varying the contributions of the 
p(770), p(1450), and p(1700) to match the data below x = 0.925, then extrapolating into 
the region 0.925 < a; < 1 sensitive to rrii,^. Q The resulting likelihood function (Fig. ^) is 
sensitive mainly to the parameters describing the p(1700). If we vary the amplitude within 
the experimental accuracy of our matching and the p(1700) mass and width within their 
experimental errors, we can raise the 95% likelihood limit by 4 MeV. 

The effects of other systematic uncertainties have been similarly evaluated. The resulting 
shifts of the 95% limit on rriy^ are listed in Table |T|. Assuming that the various effects are 
independent, we combine the limit shifts in quadrature to get a net systematic shift of 6.4 
MeV. Following the practice in reports of previous experiments on m^^ limits, we add this 
shift linearly to the raw limit from Fig. ^ 

m^^ < 28 MeV (95% confidence). 

An important difference between this measurement and u^. mass limits from previous 
e~^e~ experiments is the size of the event sample used. We observe 543 events in the sensitive 



^A Poisson coefficient expressing the dependence of the number of observed events on m^^ (as 
used in ref. |jl^) is not used here because the large number of events in the fit region of this work 
is insensitive to the neutrino mass scale in question. 

^Note that if we were to force the spectral function to agree with the data all the way to the 
kinematic end point {x = 1), we would get a (biased) 95% limit of 17 MeV for m^^. This is 
significant in that it represents the lowest limit one could a priori expect to obtain with the given 
statistical accuracy and background, assuming m^^ = 0. 

9 



portion of the kineinatically allowed region, x > 0.925. This and the fact that we include 
background in the fit lead to a limit that has little sensitivity to chance fluctuations in 
the population of individual events near the endpoint. Although this experiment has its 
own statistical and systematic uncertainties that prevent a significant improvement in the 
limit value, the analysis is quite complementary to previous low statistics experiments and 
confirms their conclusions. 



Both this measurement and the previous CLEO measurement [|14| using r — * Sttz/,- show 
a broad likelihood maximum near the higher end of the allowed range of the tau neutrino 
mass. Since such a behavior is not unlikely even if m,^^ = (as verified by Monte Carlo 
experiments), we do not regard it as significant. It does imply, however, that combining the 
results of the two CLEO measurements to make a joint likelihood curve does not significantly 
improve the mass limit. 
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TABLE I. 95% confidence upper limits on rriu^ obtained from kinematics of r decays. 



Experiment 



Ref. 



Decay 



MeV 



ALEPH 

OPAL 

ARGUS 
CLEO 



ref. flS 



ref. PI 



ref. I? I 

ref. Il4| 

this analysis 



37r± 
57r±(7rO) 
57r± 
37r± 
57r± 
57r±,37r±27r° 



25.7 

23.1 

43.2 

35.3 

31 

30 

28 



TABLE IL Fraction of the r — > /i h^h -k^u-j- data contributed by each of the signal modes, 
based on known branching fractions S. 



Decay mode 



of signal 



K-K+T:-i:^Ur 



92 
3 

1.6 
1.5 
1.1 
0.5 
0.2 



TABLE in. Systematic uncertainty sources and the shifts they induce in the rrijy^ upper limit. 



Source 



MeV 



TT^ energy scale 

Track momentum scale 

Spectral function 

qq and r background corrections 

Monte Carlo statistics 

Resolution function 



3.7 
3.3 

4.0 
0.8 
0.5 
0.4 



Quadrature sum 



6.4 



12 



> 



(0 

C 

> 

LU 



4000 



3000 



2000 



1000 



1 \ \ r 



J \ L 



1 \ \ r 







J \ L 



0050399-002 

1 \ \ \ 




J \ L 



0.5 



1.0 




M(7rV"7r°)(GeV) 

FIG. 1. Invariant mass distribution for tt^tt^tt^ combinations in data, with two entries per 
event. No backgrounds have been subtracted. 
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FIG. 2. Invariant mass of tttt combinations in events outside the lu peak: 7r"'"7r'^ (a, one entry 
per event), 7r~7r'^ (b, two entries per event), tt^tt" (c, two entries per event). The data distribution 
is represented by the points with error bars and Monte Carlo by the histogram. 
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FIG. 3. Four-pion invariant mass for events within the uj peak (a) and outside the to peak 
(b). The data distribution is represented by the points with error bars and Monte Carlo by the 
histogram. 
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FIG. 4. Distribution in scaled hadronic energy versus scaled hadronic mass for events in the 
final data sample, shown for the full range (a) and for a restricted range near the kinematic end 
point (b). The curves show the boundaries of the allowed kinematic region for r -^ TT~TT~^n~TT^Ur 
assuming niu^ = in (a) and (b) and m^^ = 30 MeV in (b). The dashed lines in (a) show the 
boundary of the fit region. 
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FIG. 5. The measured likelihood evaluated for a sequence of mi^^ values. The curve is a cubic 
spline fit through the points. The shaded area is 5% of the integral under the curve. 
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